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1 Introduction 

A bi-Hamiltonian manifold M is a smooth manifold endowed with two compatible bi-vec- 
tors P, P' such that 

[P,P] = [P,P'} = [P',P'] = 0, 

where [•, •] is the Schouten bracket. The bi- vectors P, P' determine a pair of compatible Poisson 
brackets on M, for instance 

{f(z),g(z)} = (df,Pd 9 ) ^ P ij (*) %lf% { *\ (1-1) 

and similar brackets {-,•}' to P' . 

Dynamical systems on M having enough functionally independent integrals of motion H%, . . ., 
H n in involution with respect to both Poisson brackets 

{H l ,H j } = {H i ,H j } 1 = 0. (1.2) 

will be called bi-integrable systems. 

A sufficient condition in order that integrals of motion H\, . . . ,H n be in bi-involution is 
that the corresponding vector fields are bi-Hamiltonian vector fields \13\ [T7] , which form 
a so-called anchored Lenard-Magri sequence 

PdHx = 0, X H , = PdHi = P'dHi-i, P'dH n = 0. (1.3) 

The class of manifolds we will consider are particular bi-Hamiltonian manifolds, to be termed 
o;A^-manifolds, where one of the two Poisson bi- vectors is nondegenerate (say P) and thus defines 
a symplectic form u = P~ l and, together with the other one, a recursion operator [31 0] 

N = P'P- 1 , (1.4) 



*This paper is a contribution to the Vadim Kuznetsov Memorial Issue "Integrable Systems and Related Topics" 
The full collection is available at |http://www.emis.de/journals/SIGMA/kuznetsov.html| 



Yu.A. Grigoryev and A.V. Tsiganov 



and its dual N* = P~ l P' . Operator N is called a Nijenhuis operator |131 [T7] or hereditary 
operator [TQJ[8]. 

One of the main property of the recursion operator ./V is that bi-vectors P = P^ and 
P 1 = belong to a whole family of compatible Poisson tensors 

p(k) = N k p ^ £ Lg j 

which define a family of compatible Poisson brackets {•, ■} on the wiV-manifold M. 

Another useful property of ./V is that normalized traces of the powers of N are integrals of 
motion satisfying Lenard-Magri recurrent relations (jl.3p |17j : 

1 



H j = —trN j . (1.6) 



The class of coordinates, called Darboux-Nijenhuis coordinates, are canonical with respect to u 
and diagonalize recursion operator N. According to [31 H], the n-tuple (Hq, . . . , H n ) of Hamil- 
tonians on M (where n = \ dimM) is separable in Darboux-Nijenhuis coordinates if and only 
if they are in involution with respect to both Poisson brackets (jl.2p . 

In this paper we compare two known families of the separated variables for the open Toda 
lattice [TBI 03] with the corresponding Darboux-Nijenhuis coordinates. 



2 The separation of variables method 



In the separation of variables method variables we are looking for complete integral S(q, t, ati, 
a n ) of the Hamilton- J acobi equation 



OS 



— +H[q,—,t 



Of 



dS 



dq 



0. 



det 



8 2 S 



dqidaj 



where q = (qi, . . . , q n ), in the additive form 

n 

S(q,t,ai,...,a n ) = -Ht + y]Sj(qi,ai, ■ ■ ■ ,a n )- 



(2.1) 



(2.2) 



i=l 



Here the i-th component 5, depends only on the i-th coordinate qi and n parameters a\, . . . , a n 
which are the values of integrals of motion. In such a case H is said to be separable and 
coordinates q are said to be separated coordinates for H, in order to stress that the possibility 
to find an additive complete integral of (|2.2p depends on the choice of the coordinates. 

For any complete integral of the equation (|2.ip solutions qi = qi(t,a,[3) and pi = pi(t,a,(3) 
of the Hamilton equations of motion are found from the Jacobi equations 



R. - - — 

HI — <1 ' 
OOi 



dS_ 

dqi 1 



1, .. . , ra- 



in the separation of variables method the each second Jacobi equation 



Pi = ■S-S i (q i ,a 1 ,...,a n ), 
dqi 



(2.3) 



(2.4) 



contains the pair of the Darboux coordinates pi and qi only. These equations and their more 
symmetric form 



$i(qi,Pi,ai, ...,a n ) = 0, 



det 



(2.5) 



are called the separated equations. 



On the Darboux-Nijenhuis Variables for the Open Toda Lattice 



3 



The separated variables (p, q) are defined up to canonical transformations p, — > fi(pi, qi) and 
Qi 9i{Pii Qi) and integrals of motion Hj = oej can be always replaced with Hj = (f>j{H\, . . . , H n ). 
Such transformations change the form of the separated equations (|2.5p . 

For example, we can use canonical transformation to the action-angle variables 

n 

Here W = ^ Si(qi,a\, . . . ,a n ) is Hamilton's characteristic function. The action-angle variab- 

i=i 

les /, w are the special separated variables which allow us to linearize equations of motion 

h = °> ™j = F j{h, ■■■Jn), j = 1, • • • , n. 

However, the action-angle variables are not always convenient, for instance for the quantum 
integrable system. 

Example 1. Consider a two-particle open Toda chain with the following integrals of motion 

2 2 

H 1 =p 1+P2 , fl- 2 = a+£2+ e «i-*.. (2.7) 
Variables 

«i,2 = ± 92), u>i j2 = -f=(pi ± P2) (2.8) 

are separated variables, because substituting these variables in the definition of integrals of 
motion we obtain desired separated equations (|2.5p 



2 tt2 

*i(ui,«i,a) = v^^! - Hi = 0, $ 2 (w 2 ,u 2 ,a) = e^ 2 + -| - # 2 + = 0. 
(95- 

Provided Hj = a^ and ui = -= — , these separated equations 

dvi 

2 ~ ™2 



can easily be integrated by quadratures 

<? - aiVl O 



S 2 = =F V 4a 2 - af - 4e^ 2 ± W 4a 2 - af arctanh 



'4a 2 - a\ - 4e^ v 
y 7 4a 2 — a\ 



Then we have to substitute the corresponding Hamilton principal function S = —a 2 t + S\ + S 2 
in the Jacobi equations 



01 = — — — = 7= H 1 — arctanh 



Ua 2 -a\- 4e^ 2 \ 



— — = -j= -| — arctann — 

<9ai V2 V4a 2 - I V 4a 2 - a\ I 

OS 2 ( \l±a 2 -aj-4e^ 

02 = — = t — arctanh 



da 2 A /4a 2 - a\ \ a/4o 2 



a; 
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and solve the obtained equations with respect to 

Vl = h V2y9l, 

V2 



Using these solutions and canonical transformation (12. 8p we obtain closed equations for trajec- 
tories of motion in the original variables 

qi2 = a Jf-M ± l ln ^ _M), ln L 8h (y^f « - ft)] ] - A . (2.10) 



In similar way we can get pi 2(i) from the second Jacobi equations (|2.4p . 

Parameters /Ji^ may be excluded by the shifts t — > t + 02 arid 9i — > 9i + 0\ and, thus, the 
solutions 91,2 (*) depend on two values of integrals of motion ai,2 only. 

Now from (|2.6|) one can easily find the action variables 



I 1 = a l = V2u 1 , I 2 = 4a 2 - a\ = 2u\ + 4e^ 2 (2.11) 
and the angle variables 



U2 



arctanh , , 

ds V! ds i Vv«i+ 2e ^V , 010 , 

^ = m = 7S' W ^mr-W2 jr^z • (2 - 12) 



2 

In these coordinates symplectic form has a standard form 

2 2 

uj = ^2 d-Pj A dqj = dlj A dwj 
3=1 3=1 

and equations of motion are linearized 

• . dH h . dH 1 

/j = ' W1 = ^ = T' " 2 = 572 = 4- 

Using solutions of these equations and inverted canonical transformation I,w — > p, g one can 
easily derive equations for the trajectories of motion (|2.10p . 

Separated variables V12 and u\2 can also be used in quantum mechanics, for example, to 
find the spectrum of Hamiltonian of a periodical two-particle Toda chain [1] 



2.1 Darboux— Nijenhuis coordinates 

In this section we will describe a class of canonical coordinates on cjA^-manifolds, called Darboux- 
Nijenhuis coordinates. They will play the important role of variables of separation for (suitable) 
systems on wA^-manifolds. 

By definition a set of local coordinates (xi,yi) on an wiV-manifold is called a set of Darboux- 
Nijenhuis coordinates if they are canonical with respect to the symplectic form 

n 

to = P^ 1 = dyi A dxi 
i=i 
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and put the recursion operator TV in diagonal form, 

This means that the only nonzero Poisson brackets are 

{xi,yj} = 8ij, {xi,Vj}' = XiSij. (2.14) 

The distinguishing property of the pairs of Darboux-Nijenhuis coordinates (xi,yi) is that their 
differentials span an eigenspace of N* , that is, satisfy the equations 

N*dxi = Xidxi, N*dy { = Xidyi (2.15) 

As a consequence of the compatibility between P and P' , the Nijenhuis torsion of N vanishes 

T N (X, Y) = [NX, NY] - N([NX, Y] + [X, NY] - N[X, Y]) = 0, (2.16) 

here X, Y are arbitrary vector fields on M. 

According to the Frolicher-Nijenhuis theory [9], condition ([2.16]) implies that the distribution 
of the eigenvectors of N is integrable. In application to Darboux-Nijenhuis coordinates it 
means that it is possible to find by quadratures 2n functions (xi,yi) directly from the equa- 
tions (]2.15p [111116] . Hence, the Frolicher-Nijenhuis theory allows one to construct the Darboux- 
Nijenhuis variables as solutions of the equations (|2.15p or their equivalent (I2.14|) . 

However, for any function f(xi,yi) equation 

N*df(xi,yi) = \idf(xi,yi) 

is also satisfied. Hence, any pair of Darboux-Nijenhuis coordinates is defined up to arbitrary 
canonical transformations. So in the general case equations (|2.15p have infinitely many solutions 
and, thus, it is impossible to create a common effective algorithm to solve these equations. 

This problem can be partially solved. As a consequence of the vanishing of the Nijenhuis 
torsion of N the eigenvalues Aj always satisfy (12.15P 

N*d\i = XidXi. 

The eigenvalues Aj and their conjugated variables [i{ are the special Darboux-Nijenhuis coordi- 
nates [4], which are distinguished because \ are simply the roots of the minimal characteristic 
polynomial of N 

n 

Ajv(A) = (det(7V - AI)) 1/2 = A n - (ciA™" 1 + • • • + c n ) = J](A - A,)- (2-17) 

i=i 

Here Cj(X\, . . . , A n ) are elementary symmetric polynomials of power j which are related with 
integrals of motion Hj (ll.6p by the Newton formulas. 

The complimentary variables \ij must be calculated as solutions of the overdetermined system 
of partial differential equations ([2. 15[) 

N*dfjbi = XidpLi, {Xi,fj,j} = 5ij, {m, fj,j} = 0, (2.18) 

whose solutions /ij are still determined up to infinitely many canonical transformations fij — > 
to + f( x j)- 

An effective algorithm for calculating Darboux-Nijenhuis variables has not been developed 
yet despite these variables play an important role in the method of separation of variables due 
to the following theorem. 



1. 
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Theorem 1 ([31 0]). If M is a In- dimensional uN-manifold such that in the neighborhood of 
any point z G M recursion operator N has n different functionally independent eigenvalues and 
if {Hi, . . . ,H n } is a family of independent functions on M, then the following statements are 
equivalent: 

1) functions {Hi, . . . , H n } are in bi-involution (|1.2jl : 

2) the Lagrangian foliation T defined by {Hi, . . . , H n } is separable in Darboux-Nijenhuis 
coordinates; 

3) the distribution T> tangent to the foliation defined by {Hi, . . . , H n } is Lagrangian with 
respect to uj = P^ 1 and invariant with respect to N . 

It easy to prove that for integrable after Liouville bi-Hamiltonian systems with the integrals 

of motion Hj (jl.6p recursion operator N has exactly n functionally independent eigenvalues, 

n 

because Hj = j^ 1 H an d 

dHi A • • • A dH n = JJ(Aj - Xj)dXi A • • • A d\ n . 

For bi-Hamiltonian systems integrals of motion Hi (jl.6p define a distinguished bi-Lagrangian 
foliation, called principal foliation [3]. 

For generic bi-integrable systems invariance of the distribution T> with respect to N means 
that there is a control matrix F with eigenvalues (Ai, . . . , X n ) such that 



N*dHi /•',////,. 



1. 



n. 



(2.19) 



fc=i 



In the Darboux-Nijenhuis coordinates equations (|2.19)) may be considered as the known Levi- 
Civita criterion for separability [15], see [3]. 

In particular case for bi-Hamiltonian systems (II. 3|) the control matrix F has the following 
form 



/o 





1 







1 



■■■ 0\ 
... o 

1 

cij 



(2.20) 



Here c/% are coefficients of the characteristic polynomial Aat(A) (|2.17p of the recursion opera- 
tor N . 

Example 2. According to [21 [6], for a two-particle open Toda chain recursion operator reads as 



/ 


pi 








1 \ 







P2 


-1 










_ e <?i-92 


Pi 





V 


g<?l-<?2 








P2 J 



N 



Integrals of motion (|2.7D are reproduced by f 1 1 . 6 [> . Therefore eigenvalues of N 
P1+P2 . \Jp\+p\- 1p\P2 + 4exp(gi - q 2 ) 



+ 



are a half of the special Darboux-Nijenhuis coordinates for the open Toda lattice. Of course, 
these coordinates are the action variables A,- = 0. 
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The complimentary variables fi± t2 must be calculated as solutions of the overdetermined 
system of partial differential equations (|2.18|) . However even in this simple case we could not 
directly solve these thirteen PDEs (|2.18p for the two unknown functions H-i(qi,q 2 ,Pi,p 2 )- 

Nevertheless we can find variables fj,i by using the action-angle variables p. lip . (12.121) ob- 
tained before. Since 

A, = ^4^, (2.21) 
the second half of the special Darboux-Nijenhuis coordinates reads as 

2 ^T 2W2 = T arctanh Pl ~ P2 . (2.22) 

2 VvG°i -P2) 2 + 4exp(gi - q 2 ) ) 

By definition action variables Aj are roots of the minimal characteristic polynomial 

Ajv(A) = A 2 - (pi + p 2 )\ + Pl p 2 - e« 
whereas angle variables fii can be defined in the following way: 

fM = lnB(\ i ), B(\) = -e^(X- Pl ). 

Summing up, to build the special Darboux-Nijenhuis variables we de-facto had to use another 
set of separated variables Vi, Ui which are not Darboux-Nijenhuis variables. 

In the next section we discuss how analogous variables u and v are related with the special 
Darboux-Nijenhuis variables for n-particle open Toda lattice. 

3 Open Toda lattice 

Let us consider open Toda associated with the root system of srf n type. The Hamilton function 
is equal to 



1 n n—1 



i=l i=l 

Here p, q are Darboux coordinates on the manifold M ~ M? n 

{q i ,p j } = S ij , {pi,Pj} = Ui,qj} = 0. (3.1) 

Bi-Hamiltonian structure of Toda chains was investigated both in terms of physical variables 
(p, q) and in terms of so-called Flaschka variables [7] . We will use original physical variables p 
and q, for which second Poisson tensor has the form [21 [6] 

Tt 1 r\ r\ TL r\ r\ Th r\ r\ 

p> y 1 9i-9i+i A ^ A ^ A (3 2) 

dp i+1 dpi ^ Pt dqi dpi dqj dq{ 

Throughout the rest of the paper N is recursion operator N = P'P^ 1 for the open Toda 
lattice, where P' is given by (|3.2p and P is canonical tensor associated with canonical Poisson 
brackets (ED in R 2n . 
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3.1 The Moser method 

In this section we will give a brief review of a method for integrating equations of motion for an 
open Toda lattice proposed by Moser [18]. A contemporary algebraic-geometric review of the 
Moser method can be found in [25]. 

Let us start with a n x n Lax matrix L with the components 

Ljk = Pj5j,k + e 2 (Sj tk+1 + £j-i,fc). (3.3) 

Following [18] we define the Weyl function 

X(X) = (R{X)a, a), a = (0, 0, . . . , e qn/2 ), (3.4) 

where R(X) = (L — A/) -1 is a resolvent of the Jacobi matrix (|3.3p and a is dynamical normaliza- 
tion of the corresponding Baker-Akhiezer function. It is known since the work of Stiltjies [22], 
the Weyl function plays the key role in reconstruction of the matrix L from its spectral data 
and it is a ratio of the two monic polynomials 

*w = fgf (3.5) 

where B(X) is a polynomial of degree (n — 1) and A{X) = det(L — XI) is a polynomial of degree n 
with distinct roots Xj 

n 

A(X) = Y[(X-Xj), B(X) = -e 9 "(A n ~ 1 + b 2 X n ~ 2 -\ \-b ). (3.6) 

i=l 

Proposition 1. For the open Toda lattice two Poisson brackets associated with tensors P = P^ 
(jl.5|) and P' = P^ ()3.2[) form quadratic algebras for the polynomials A(X), B(X) and for the 
Weyl function A'(A). 

It is easy to prove that at k = 0, 1 these quadratic brackets are equal to 

{A(X),A(^ ={B(X),B(»)}W =0 , (3.7) 

{A(X),B(^ 



H k A(X)B(fi) - X k A(n)B(X) 



X — n 
and 



X — fi 



(3. 



At k = brackets (|3,7p are the part of the Sklyanin brackets [19] and nontrivial expression in 
the right hand side of (|3.7p is called a Bezoutian of polynomials A and B [12]. 

At k = brackets (|3.8|) give the Atiyah-Hitchin Poisson structure in the space of meromorphic 
maps X(X) : CP 1 -> CP 1 [l]. 

Remark 1. According to [5] the Poisson brackets between polynomials -A(A) and B((i) associa- 
ted with Poisson tensors P^ (|1.5p at A; = 0, 1, . . . , n are equal to 

{m Bm m = i(A)#iM-iMB"i(» 

A — /x 

= n k A(X)BQi) - X k A^)B{X) + ^^(^W^M^ 
X — fj, X — fj, 
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The right hand side of (|3.9p contains either remainder 

B^(X) = X k B(X) modA(A), (3.10) 
either result of simple division in the space of polynomials 

*« = ^, (3-n) 

which is polynomial part of the Laurent decomposition of the ratio X k B(X)/A(X). It is easy to 
see that /?W = and 0& = -e q " and from fl3j| one gets fl37 



In the Moser approach to the open Toda lattice [18] we can introduce the action-angle vari- 
ables using the Weyl function (|3.4p . Namely, action coordinates are poles Aj (|3.6p of the Weyl 
function A? (A), whereas angle variables are given by 

fij = In B(Xj). 

It is easy to check that polynomial A(X) = det(L — XI) = Ajv(A) is a minimal characteristic 
polynomial of the recursion operator N. Moreover, it follows from (|3.7p that Aj and fij satisfy 
the necessary relations (I2.14P 

We could avoid calculating of the second Poisson brackets between A and B by means of the 
equation 

N*dA N (X) = XdA N (X) + A N (X)dd, (3.12) 

which can be obtained from the recurrent Lenard-Magri relations (jl.3p rewritten in the form 
N*dHk = dH^ + i and from the Newton formulas connecting with coefficients [3]. Com- 
bining this equation with a result of [23] 

B (A) = - e ^A Ajv (A) (3.13) 
aci 

one can easily prove that B{X) is a Stackel function [JJ 

N*dB(X) = XdB(X) + A 7V (A)de 9 ™ 
and, thus, 

N*dfj,i = Xidm. 

Summing up, the Moser variables are the special Darboux-Nijenhuis variables for the open 
Toda lattice. The corresponding equations of motion have the form 

{Hi, Xj} = d n Xj = 0, {Hi, fij} = d Ti fij = A*" 1 . 

Evolution of the variables fij with respect to the times Tj conjugated to the bi-Hamiltonian 
integrals of motion Hi ()1.6p is linear. 

Brackets (|3.8p for the Weyl function are invariant with respect to linear- fractional transfor- 
mations 

and, therefore, we can introduce another family of the separated variables using the same Weyl 
function [2S]. This second coordinate system is considered in the next section in framework of 
the classical r-matrix theory. 
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3.2 The Sklyanin method 

In this section we briefly discuss an application of the generic Sklyanin method |X9|, [20] to the 
open Toda lattice. Let us consider a 2 x 2 monodromy matrix 



T(A) 



A B 
C D 



(A) = Li(A) • • • L n -i(X)L n (X), 



(3.15) 



where 



u 



A -Pi 



-e qi 




Monodromy matrix T(A) (|3.15|) is the Lax matrix for periodic Toda lattice [TU], whereas the 
Lax matrix for open Toda lattice is equal to 



T (A) = KT(X) 



A B 




(A), 



K 



1 




The entries A and B of this 2x2 Lax matrix coincide with denominator and numerator of the 
Weyl function ()3.5[) respectively. 

Proposition 2. The Poisson brackets between entries of the monodromy matrix T(A) (|3.15|) 
associated with Poisson tensors P( fe ) (jTj5j) at k = 0,1 are equal to 

{r(A),T(M) } W = r 12 (A, /x)r(A)T(M) -T(A)T(/x) r 21 (A, M ) 
+T(A) s(\, m )T(m) -T(A) s(A, /i)T(^), 
where first r-matrix is equal to 



(3.16) 





/ 


1 








\ 


-1 







1 A fc +M fc 
1 2 


// 





A — 







A fc 


1 2 







I 











1 / 



whereas the second matrix s(A, /x) reads as 
s(A,m) 





/ o 











\ 


-1 





\ k -tl k 

2 










-A* 








\ k -fl k 

2 









1 











/ 




2 











Here, T(X) = T(X) /, T(/i) = I T(/x) and II is a permutation matrix in auxiliary space, i.e. 
HX Y = Y All /or arbitrary matrices X, Y 

At A; = first matrix r(A,/i) = —(A — /i) _1 II is the standard rational r-matrix, while the 
second matrix s(A,//) is equal to zero. In this case quadratic brackets (|3.16p are the Sklyanin 
bracket OS]. 



Proposition 3. In generic case the Poisson brackets {-,-}^ associated with Poisson tensors 
?W (dSJ at k = 0,1,... , n are equal to: 

{A(\),A(ti)}W = {B(X), Bfa)}M = {C(X),C(^ = {D(X), =0, 
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X — fi X — fi 

X — fi X — fi 

{fl (A), BM} <») = yWyW) _ ^-^(AMM, 

A — // 

/j k A(X)D(fi) - X k A(fi)D(X) 



{B(X),C(fi)}^ 

and 

{A(X),D^)] 



X — fi 

(k) _ X k {C(X)B(fi)-C(fi)B(X)) 
X — fi 

B{fi) (pW (A) - /3 W - ( 7 [fc] (A) - 7 [fc] (/*))). 



A — /U 

X — fi 



Here fffl = X k B(X)/A(X) (|3.1ip and = A fe C(A)/ J 4(A) are polynomial parts of the Laurent 
decompositions of quotients of the corresponding polynomials. 

So, at k > 1 we have to add to the matrices ri 2 and s dynamical terms proportional to the 
functions 0^ and 7^. 

According to [7[ [191 HO] the separated coordinates are poles of the corresponding Baker- 
Akhiezer function with the standard normalization a = (0, 1). In this case the first half of 
variables are coming from (n — 1) finite roots and logarithm of leading coefficient of the non- 
diagonal entry of the monodromy matrix 

n-i 

B(X) = -e~ Vn H(X- Uj ), v n = -q n , (3.17) 
Another half is given by 

n 

Vj = —\uA{uj), j = 1, . . . , n — 1, and u n = — ci = ^^Pi- (3.18) 

i=l 

In the separated variables polynomial A(X) reads as 

(n \ n— 1 n— 1 n— 1 , 

a + e «i n ( a - «*) + e ^ n ( 3 - 19 ) 

This definition of the separated variables is obviously related with the following transformation 
X -» Af" 1 of the Weyl function [25]. 
It follows from pTTBj) that at fc = 0, 1 

{.4(A), 5(^)}( fc ) = ^— (fi k A(X)B{fi) - X k A(fi)B(X)) . (3.20) 

Substituting these brackets into the equations 



= {A(\) i B(u j )}M = {A(X),B(fi)}^ +B'( Uj ){A(X), Uj } 



(k) 
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one gets 

n— 1 , 

{A(\), Uj }W = \ k A(uj) J] — ^, j = 1, . . . ,n - 1. (3.21) 

. , . Uj Ui 

If A = Ui this implies part of the necessary relation (|2.14p 

{uj,Vi}(V =u^Sij, k = 0,l, i,j = l,...,n- 1. 
Now we can collect all the coefficients with the same powers of A and u in (|3.2(jp and prove that 

{u n ,v n } = l, {ui,u n } = 0, {v n ,Vi} = 0, 

{u n ,V n }' = - 2jlii, {u n ,Ui}' = fi = — , {v n ,Vi}' = -1. 

In a similar manner from {B(\),B(n)}W = one gets 

{ui,Uj}( k) = {v n ,Uj}W = 0, i,j = 1, ...,n-l, 
and from {A(A), A(/i)}( fc) = it follows that 

{ii n ,Ui} = 0, = 0, i,j = 1, . . . ,n - 1, 



{u n ,Vi}' = gi = -e 



B'(X) 



{vi,VjY = 0. 



X=ui 



Here ^4'(A) and B'(X) are derivatives by A. 

Summing up, the 2n separated variables V{ and Ui are the Darboux variables 



uj = 

i=l 



dttj A (if j , 



but the corresponding recursion operator N consists of two diagonal and three non-diagonal 
terms 



( d d \ f d d 

N = \^ Ui ( — ® dUi + — (g> dVi I - Ui I - — <g) cfe n + - — ® d« n 
f-f V^ u i ^ / ■ , \ou n dv n 

1=1 1=1 

n-l 



+ ^ ffi ( 7T- ® d«n + j— ® dVi \ . (3.22) 



i=l 



Thus 2n separated variables Vi and Uj obtained in framework of the Sklyanin method [19] are not 
the Darboux-Nijenhuis variables. Evidently it is related with the difference in the form of the 
corresponding separated equations, which follow directly from the definitions of Vj and u n (13. 181) 

$j = e~ Vj — A(uj, a\,..., a n ) = 0, = u n + ax = 0. (3.23) 
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The first (n—1) equations of motion, see (|3.2ip . are linearized by the Abel transformation [19^ 121) 
f,,,,^P 1 X7 _i X^dX 



k=l 



Atv(A) 



where {crj} is a basis of Abelian differentials of first order on an algebraic curve z = Ajv(A) 
corresponding to separated equations (j3.23[) . 

It means that we can introduce the action-angle variables 



Ij = Cj, w\ = v n , and w j+i = 



ra-1 rVk 



k=l 

such that evolution of variables Wj with respect to times Tj conjugated to the action variables Ij 
are linear. These action variables are related with the previous ones 

n n 

A(X) = H(X-X i ) = X n -Y,Ij>> n - j - 

i=l j=l 

Bi-Hamiltonian integrals of motion Hj ()1 .6|) produce the flows which preserve the spectrum of 
the n x n Jacobi matrix L (13. 3| . Integrals of motion Ij = Cj produce the transversal to the 
isospectral manifolds flows, which preserve the divisor [25j . 

As sequence the special Darboux-Nijenhuis variables A, [i are dual to the Sklyanin variab- 
les u, v. Namely, Aj, m are roots of polynomial ^4(A) and values of polynomial B(X) at A = Xi, 
while Uj, Vj are roots of polynomial B(X) and values of polynomial A(X) at A 



j- 



Remark 2. From the factorization of the monodromy matrix T(A) tj3. 15 j> one gets 

B n {X) = -e«»4v-i(A) => B n (uj) = -e^An^Xj) = 0. 

This implies that the Moser variables Xj for a (n — l)-particle chain, i.e. special Darboux- 
Nijenhuis variables, coincide with the Sklyanin variables Uj, i = 1, ...,n— 1 for a n-particle 
chain. 

We can prove this fact directly using matrix representation of the recursion operator for 
(n — l)-particle chain, which can be obtained by the matrix representation (I3.22p deleting the 
n-th and 2n-th rows and columns. 

Remark 3. The Sklyanin variables are the separated variables for open and periodic Toda 
lattices simultaneously, in contrast with the Moser variables. For the periodic Toda lattice the 
Darboux-Nijenhuis variables were constructed by using Flaschka variables in [3]. 

Example 3. At n = 2 the Sklyanin variables are equal to 

vi = -ln(-e qi ~ q2 ), v 2 = -q2, m = p 1 , u 2 = pi + jp 2 , 

and 



N 



2ui 











— e 




- u 2 


U 2 - Ui 


g-«l 










-1 


Ul 


2«i - 


- u 2 


1 








u 2 - 





We have three different families of the separated variables for the Toda lattice at n = 2 only. 
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4 Conclusion 

For the open Toda lattice associated with the root system of s$ n type we prove that the Moser 
variables are special Darboux-Nijenhuis variables, while the Sklyanin variables are "almost 
Darboux-Nijenhuis variables", in which recursion operator consists of a diagonal part and two 
non-diagonal rows and columns only. 

The similar results for the generalized open Toda lattices associated with the root systems 

of 

and Sl n type were found in [24j . 
Nevertheless we have to underline that construction of the Moser variables and the Sklyanin 
variables directly in the framework of the bi-Hamiltonian geometry is an open question. 

Acknowledgements 

We would like to thank I.V. Komarov and V.I. Inozemtsev for useful and interesting discussions. 
The research was partially supported by the RFBR grant 06-01-00140. 

[1] Atiyah M., Hitchin N., The geometry and dynamics of magnetic monopoles, M.B. Porter Lectures, Princeton 
Univ. Press, 1988. 

[2] Das A., Okubo S., A systematic study of the Toda lattice, Ann. Phys., 1989, V.30, 215-232. 

[3] Falqui G., Magri F., Pedroni M., Bihamiltonian geometry and separation of variables for Toda lattices, 
J. Nonlinear Math. Phys., 2001, V.8, 118-127, jnhn.SI/0002008 

[4] Falqui G., Pedroni M., Separation of variables for bi-Hamiltonian systems, Math. Phys. Anal. Geom., 2003, 
V.6, 139-179, [rilm.SI/0204029| 

[5] Faybusovich L., Gekhtman M., Poisson brackets on rational functions and multi Hamiltonian structures for 
integrable lattices, Phys. Lett. A, 2000, V.272, 236-244, [m"in.SI/0006045| 

[6] Fernandes R.L., On the master symmetries and bi-Hamiltonian structure of the Toda lattice, J. Phys. A: 
Math. Gen., 1993, V.26, 3797-3803. 

[7] Flaschka H., McLaughlin D.W., Canonically conjugate variables for the Korteweg-de Vries equation and 
the Toda lattice with periodic boundary conditions, Progr. Theoret. Phys., 1976, V.55, 438-456. 

[8] Fokas A., Fuchssteiner B., Symplectic structures, Backlund transformations and hereditary symmetries, 
Phys. D, 1981, V.4, 47-66. 

[9] Frolicher A., Nijenhuis A., Theory of vector- valued differential forms. I. Derivations of the graded ring of 
differential forms, Nederl. Akad. Wetensch. Proc. Ser. A, 1956, V.59, 338-359. 

[10] Gelfand I.M., Dorfman I.Ya., Hamiltonian operators and algebraic structures associated with them, Funkt- 
sional. Anal, i Prilozhen., 1979, V.13, N 4, 13-30. 

[11] Gelfand I.M., Zakharevich I., On the local geometry of a bi-Hamiltonian structure, in The Gelfand Mathe- 
matical Seminars 1990-1992, Editors L. Corwin et al., Boston, Birkhauser, 1993, 51-112. 

[12] Krein M.G., Naimark M.A., The method of symmetric and Hermitian forms in the theory of the separation 
of the roots of algebraic equations, Linear and Multilinear Algebra, 1981, V.10, N 4, 265-308 (translated 
from the Russian by O. Boshko and J.L. Howland). 

[13] Kosmann-Schwarzbach Y., Magri F., Poisson-Nijenhuis structures, Ann. Inst. Poincare (Phys. Theor.), 
1990, V.53, 35-81. 

[14] Komarov I.V., Tsiganov A.V., Two particle quantum periodic Toda lattice, Vestnik Leningrad Univ., 1988, 
V.2, 69-71. 

[15] Levi-Civita T., Integrazione delle equazione di Hamilton- Jacobi per separazione di variabili, Math. Ann., 
1904, V.24, 383-397. 

[16] Magri F., Geometry and soliton equations, in La Mecanique Analytique de Lagrange et son heritage, Atti 
Acc. Sci. Torino Suppl, 1990, V.24, 181-209. 

[17] Magri F., Eight lectures on integrable systems, in Integrability of Nonlinear Systems (1996, Pondicherry) , 
Lecture Notes in Phys., Vol. 495, Berlin, Springer, 1997, 256-296. 

[18] Moser J., Finitely many mass points on the line under the influence of an exponential potential - an 
integrable system, in Dynamical Systems, Theory and Applications (1974, Rencontres, BattelleRes. Inst., 
Seattle, Wash.), Lecture Notes in Phys., Vol. 38, Berlin, Springer, 1975, 467-497. 



On the Darboux-Nijenhuis Variables for the Open Toda Lattice 



15 



[19] Sklyanin E.K., The quantum Toda chain, in Nonlinear Equations in Classical and Quantum Field Theory 
(1983/1984, Meudon/Paris), Lecture Notes in Phys., Vol. 226, Berlin, Springer, 1985, 196-293. 

[20] Sklyanin E.K., Separation of variables - new trends, in Quantum Field Theory, Integrable Models and 
Beyond (1994, Kyoto), Progr. Theoret. Phys. Suppi, 1995, V.118, 35-60, |solv-int/950400f| 

[21] Smirnov F.A., Structure of matrix elements in quantum Toda chain, J. Phys. A: Math. Gen., 1998, V.31, 
8953-8971, |math-ph/9805011| 

[22] Stieltjes T., Recherches sur les fractions continues, Ann. de Toulouse, VIII-IX, 1894-1895. 

[23] Tsiganov A.V., On the invariant separated variables, Regul. Chaotic Dyn., 2001, V.6, 307-326. 

[24] Tsiganov A.V., On the Darboux-Nijenhuis coordinates for the generalizaed open Toda lattices, Theoret. 
and Math. Phys., submitted. 

[25] Vaninsky K.L., The Atiyah-Hitchin bracket and the open Toda Lattice, J. Geom. Phys., 2003, V.46, 283- 
307, | |math-ph/0202047l 



